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Transportation of Cooper-pairs by a movable single Cooper-pair-box placed between two remote 
superconductors is shown to establish coherent coupling between them. This coupling is due to 
entanglement of the movable box with the leads and is manifested in the supression of quantum 
fluctuations of the relative phase of the order parameters of the leads. It can be probed by attaching 
a high resistance Josephson junction between the leads and measuring the current through this 
junction. The current is suppressed with increasing temperature. 



Delocalization of electrons due to tunneling between 
two bodies constituting a system will, according to the 
Heisenberg uncertainty principle, lower the total energy 
of this system and, hence, lead to interbody coupling. If 
the bodies are superconductors, charge exchange due to 
Cooper pair tunneling (Josephson tunneling) causes such 
coupling. Since the Cooper-pair exchange is sensitive to 
the relative phases of the superconducting order parame- 
ters of the bodies this coupling leads to the establishment 
of phase coherence i.e. a phase ordering of the relative 
phases occurs. The characteristic energy responsible for 
this ordering is the Josephson coupling energy Ej. 

If charge exchange is accompanied by charge accumu- 
lation Coulomb energy cannot be neglected. In nanos- 
tructures, where the excess charge can be strongly con- 
fined, this may result in forces strong enough to cause me- 
chanical deformations of the systems. Coupling between 
mechanical degrees of freedom and electronic degrees of 
freedom can then arise in systems having a rigidity com- 
parable to these Coulomb forces. 

That such electromechanical coupling can lead to me- 
chanical transportation of charge in a movable nanoclus- 
ter was suggested in Q as an unusual mechanism for 
shuttling of electric charge between two bodies (see also 
1 2, pi). This triggered a number of both experimental 
and theoretical activities including where 
the possibility of coherent transfer of Cooper pairs by a 
movable superconducting island was proposed. 

There, a non stationary quantum mechanical process 
involving a movable qubit, a coherent superposition of 
differently charged states on a movable superconducting 
grain, was shown to be responsible for a finite supercur- 
rent between two remote phase coherent superconduc- 
tors. However, when confronted with this result one is 
faced by a most fundamental question: Could such a non- 
equilibrium quantum mechanical process serve as a source 
for the creation of phase coherence (phase ordering) if the 
two superconductors were initially in states with definite 



number of particles ? This paper gives a positive answer 
to the above question. 
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FIG. 1. Schematic diagram of the system discussed in 
the text. A movable nanoscale superconducting grain which 
may accomodate excess charge carried by either or 1 extra 
Cooper pairs is placed between two finite isolated supercon- 
ducting leads. Close to the leads, electrostatic gates ensure 
that the Coulomb blockade of tunneling is lifted. After a se- 
quence of grain rotations, the grain moving periodically con- 
tacting the two leads sequentially, a probe (Josephson) junc- 
tion is connected through which a current / is measured. 

To demonstrate this, we consider the same system as 
in but with the low impedance link connecting the two 
leads removed. This is shown in Fig. |l|. Initially the two 
superconducting leads, labelled left (L) and right (R), are 
in projected BCS-states and \Nfj} with Nl and Nu 
extra Cooper-pairs respectively. In terms of BCS-states 
10) with definite phase (j) we have 
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Transportation of Cooper-pairs between the leads is then 
provided by a movable nanometer-sized grain | pO| placed 
between the two electrodes. When the grain is close to 
the superconductors electrostatic gates ensure that there 
exists two nearby charge states, differing by one Cooper 
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pair on the grain, for which the difference in charging en- 
ergy l^Ec is much smaller than the Josephson energy Ej. 
Due to the smallness of the grain, Coulomb blockade of 
tunneling will restrict the available states of the grain to 
only these two states which we label \n = 0) and |n = 1), 
i.e., a single-Cooper-pair box situation is realized close 
to each lead. In this situation, when Ej ^ ^Ec: a co- 
herent superposition, a |0) -I- /? of differently charged 
states can be created on the grain. This was recently 
demonstrated experimentally Jill . 

Consider now the situation when the system is ini- 
tially in a pure state with charge neutral components 
l^i) = |n = 0) \Nl = 0) \Nr = 0). If the grain executes 
periodic motion, performing repeated alternating con- 
tacts with the two leads. Cooper-pair exchange between 
the leads will result. Thus the initial state will evolve 
into a superposition of states with different numbers of 
particles on the leads, 

\n - 0) \Nl - 0) \Nr = 0) ^ 

n=0,l Nl Nr 

The Kroenecker delta makes sure that the total num- 
ber of particles is conserved. To investigate whether this 
corresponds to a phase ordered state with a well defined 
phase difference, A$ = $fl — between the leads a 
small probe junction, a Josephson weak link, is connected 
as shown in Fig. 1. This provides the means for measur- 
ing such order since the average current / through this 
link is related to A$, 

/ = Trp/, / = /,sin($fl-$L). 

Here p is the system density matrix, 1^ is the critical 
current of the probe junction and $L,fl are the phase 
operators of the respective lead. 

To be more concrete, consider the following Hamilto- 
nian describing the system dynamics 
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^ E^j(x{t)) cos($, -^) + AEc{x{t))h. 
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Hq depends on time through the grain position x(t). The 
first term represents the coupling between the grain and 
the leads in terms of the phase operators of the leads, 
^L,Ri and the grain, 0, and time varying Josephson cou- 
pling energies Ej'^{x{t)). The second term, contain- 
ing the operator n for the grain, accounts for the dif- 
ference in charging energy, AEc, associated with dif- 
ferent charge states when the grain is moving between 
the leads. In order for a phase ordering to occur two 
additional ingredients are required: First, in its present 
form [Hq , /] = which means that no sign of any or- 
dering can be detected. Second, a term describing dissi- 
pation/decoherence is needed in order for the system to 



reach a well defined final state irrespective of initial con- 
ditions. If one takes into account that any real isolated 
superconductor has a finite capacitance the first point 
no longer holds. To mimic the behavior when the leads 
have finite capacitance we restrict the Hilbert space of 
the leads to states \Nl,r — m) with |m| < N. To in- 
clude the effect of the environment we consider the effect 
of the coupling between the charge on the grain and the 
fluctuations of the gate voltage through the term 

M 

V = e(2en) ^ C^x^, < e < 1, 



where e is a dimensionless coupling constant |12|. The 
Xrv are the coordinates of a harmonic oscillator bath 
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characterized by a spectral density 



where R is the resistance and C the capacitance as seen 
from the gate. Hence, the Hamiltonian describing sys- 
tem, environment and interaction between the two is 

H{t) ^ Ho{t)+H£ + V. 

Writing the Liouvillc-von-Neumann equation in the in- 
teraction representation, 

V can be treated perturbatively to second order. Trac- 
ing out the environmental states yields for the reduced 
density matrix of the system, psit) = Ti£p{t), 

Ps{t) = Ps{to) 



dti[n{ti), / dt2Ki{ti-t2)[n{t2),Ps{h)]] 



to to 
f tl 



+ i / dti[fi{ti), / dt2K2{ti-t2){n{t2),ps{to)}] 



(1) 



Ps{to) + C[p{to)]{t). 



The kernels K12 appearing here are related to the spec- 
tral density and the inverse temperature {3 — {kBT)~^ 
through 



Kiit) — 2~ / duj J (Lu) COS Lut coin 
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and 



K2{t) 
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— — / dLoJ{uj) sinujt. 
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To find the (quasi) stationary state of the system eq. (|l|) is 
iterated. To be more precise, suppose the grain performs 
rotations with a period 4To. Introduce the following time 
labels 

• Iq: The grain leaves the left electrode. 

• ti. The grain arrives at the right electrode. 

• The grain leaves the right electrode. 

• ^3: The grain arrives at the left electrode. 

• ^4 = io + the grain leaves the left electrode. 

The reduced density matrix after a complete rotation is 
then found by (numerical) iteration, i.e. 

Ps{ti+l) = Psiti) + t[p{ti)]{ti+i). 

For the algorithm to converge the coupling to the bath 
must be weak enough. The relevant strength of the 
perturbation is characterized by the factor e^rj where 
7? = (4e2To)/(^nC) = 4T^/RqC. 

Inspired by Q we consider the grain oscillating in such 
a way that it contacts each of the leads for Tq = 1 ns and 
that it takes the same time for it cross the gap between 
the leads. Further, a reasonable value of the Josephson 
energy during the contacts is ||ll|] Ej — TfieV and we will 
assume that AEc is of the same order. For simplicity we 
take the time dependence of Ej and AEc to be step like. 
Thus close to the leads Ej takes a constant value while 
AEc is zero and vice versa. More relevant than the ex- 
act value of AEc are the dynamical phase differences x± 
picked up by the grain during the motion between the 
leads. More precisely we define 



dtAEc{t), X- = / dtAEc[t) 
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The dependence of the probe junction current on the 
dynamical phase difference is shown in figure |^. Here the 
maximum capacity of the leads is iV = 8 and the temper- 
ature set to 1 mK. Starting from an initially pure state 
the current stabilizes to a fixed value after a large num- 
ber of rotations. This current shows a distinct oscillatory 
dependence on the dynamical phases x± • The number of 
rotations required before the system stabilizes depends 
on the strength of the coupling to the bath. Since our 
approach is perturbative a very weak coupling has been 
used in order to obtain the numerical results. This im- 
plies in our case, as is evident from Fig. ^, that approx- 
imately 10'' rotations are neccessary for the value of the 



probe junction current to stabilize. In a real situation 
this time may be considerably shorter. 




FIG. 2. Average current I through a probe juction at- 
tached as in Fig. hi after many grain rotations. The magnitude 
of the current is shown as a function of the phases x± (in units 
of tt). Bright areas correspond to large current while + and 
— indicate its direction. 

Although a finite current will flow through the weak 
link in the probe junction it is not immediately clear that 
the state built up by the rotations is a BCS-type state. 
To see that this is indeed the case consider the BCS-type 
phase-states 
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Any state of the system can then be represented in the 
basis of the phases of the leads I^l,^) and the grain \(f)) 
but also in terms of the relative phases A$ = — 
and A0 = — 



|A$,A0) 



N-n 

E e~^^^^^e-^-^^\n)\N,)\~N^-n). 



,1=0.1 Nl = -N 



Since the relative phase A<f> between the grain and the 
right lead necessarily suffers from large quantum fluctu- 
ations we define the 'phase difference probability density ' 
/(A$) as the average 

f{A^) = j d{A^) (A$,A(/.|p|A$,A<^) 


In Fig. H /(A$) has been calculated as a function of the 
number of rotations performed for a system with iV = 20 
instead of TV = 8. In this simulation the fact that a defi- 
nite phase relation is built up is numerically verified. In 
the beginning of the simulation when the system is in a 
pure state with definite number of particles on each of 
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the individual components the phase uncertainty is max- 
imal, i.e. /(A$) = (27r)~^. As the system approaches a 
stable regime a sharp peak appears around a well defined 
value of A$. The location of the peak then defines the 
phase difference between the two leads. 
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FIG. 3. Phase difference probability density /(A$) as a 
function of the number of grain rotations. The graph shows 
how the system evolves from a state with maximum to a state 
with minimum quantum fluctuations as the number of rota- 
tions increases. 

The exact value of the built up phase difference de- 
pends on the parameters of the system. Especially it 
depends on the phases x± but also on the time spent 
by the grain in contact with the leads. In general, the 
final distribution /(A<i>) has two peaks separated by a 
distance tt. The relative magnitude of the peaks vary. 
For instance, the black lines, representing zero current, 
in Fig. H corresponding to the lines x+ + X- = constant 
are due to a superposition of two equally strong built 
up phases differing by tt and hence carrying an equal 
amount of current in opposite directions. A remnant of 
such a superposition can bee seen in the distribution in 
Fig. ||. Here, however, only one of the two states is stable 
while the other one rapidly decays. The largest magni- 
tude of the current is obtained when one of the two peaks 
is completely suppressed. We note further, that it is not 
possible to obtain the value of the built up phase dif- 
ference by equating the final expression for the current 
in 1^ to zero since that expression was derived under the 
explicit assumption that the phases were fixed. 

By varying the temperature, T, in the simulations the 
current I{T) through the probe junction was calculated. 
In Fig. ^ I{T)/Ic is shown as a function of temperature. 
The decrease in current with increasing temperature is 
due to an increased width of the phase distribution. At 
very low temperatures the current saturates due to the 
fact that quantum fluctuations dominate over the ther- 
mal fluctuations. Then, with increasing temperature, the 
current decays exponentially at first but with a crossover 
to algebraic decay for high temperatures. 
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FIG. 4. Probe junction current versus temperature. The 
probe junction current is here shown for flxed dynamical 
phases for a series of different temperatures. 

In conclusion we have shown that mechanically assisted 
transportation of Cooper pairs between initially uncou- 
pled, remote superconductors leads to a supression of 
quantum fluctuations of the relative phase of their order 
parameters. The movable single Cooper-pair-box, (which 
may be in a superposition of two charge states (qubit)) 
executing sequential tunneling contacts with each of the 
superconductors, is responsible for creating the strong 
entaglement causing coherent coupling. This coupling is 
supressed with increasing temperatures. 
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